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A connection is made between the random turns model of vicious walkers and random per- 
mutations indexed by their increasing subsequences. Consequently the scaled distribution 
of the maximum displacements in a particular asymmeteric version of the model can be 
determined to be the same as the scaled distribution of the eigenvalues at the soft edge of 
the GUE. The scaling of the distribution gives the maximum mean displacement /i after t 
time steps as fi = {2t) 1 / 2 with standard deviation proportional to /i 1 / 3 . The exponent 1/3 
is typical of a large class of two-dimensional growth problems. 

Non-intersecting (vicious) random walkers were introduced into statistical mechanics || as 
models of domain walls and wetting in two-dimensional lattice systems, and have also received 
attention as exactly solvable systems j?], ||, |9|, ||, [|. They can be viewed as directed lattice paths 
which start at sites say on the x-axis and finish on sites on the line y = n, with the additional 
constraint that the paths do not touch or overlap. Alternatively, vicious random walkers can be 
described as the stochastic evolution of particles on a one-dimensional lattice, which at each tick 
of the clock move to the left or to the right with a certain probability, subject to the constraint 
that no two particles can occupy the one lattice site. 

Our interest is in the random turns vicious walker model |], S|. Here, in the stochastic 
evolution picture, at each time step t (t = 1, 2, . . .) one particle is selected at random and moved 
one lattice site to the left with probability W—i, or one lattice site to the right with probability 
w\ (iw-i + w\ = 1). However, if the lattice site to the left (right) is already occupied, then the 
chosen walker moves to the right (left) with probability one unless this lattice site too is occupied. 
In the latter situation another walker is selected at random and the procedure repeated until 
one walker has been moved. The move of precisely one walker so determines the state at time 
interval t. An example of some typical configurations in the directed paths picture is given in 
Figure |l[ We remark that the random turns vicious walker model can also be regarded as a 
particular asymmetric exclusion process [17|. 

Two aspects of the theory of the random turns vicious walker model are the subject of this 
Letter. The first concerns the number of configurations that have the p walkers initially on 
adjacent lattice sites (/ = 1, . . . ,p) on the one-dimensional line, and have the walkers returned 
to the same sites after 2n time steps (for an odd number of time steps this is not possible - 
thus the reason for 2n). This will be shown to be simply related to the number of permutations 
of {1,2, . . . ,n} such that the maximum increasing subsequence has length no greater than p. 
Then known results concerning the distribution of the maximum increasing subsequence of a 
random permuation will be used to determine the distribution of the maximum displacement of 
the walkers starting near the I = 1 boundary of the initial positions for a variant of the random 
turns model. In this variant the number of walkers is greater than or equal to n, and all walkers 
must move to the left for time steps up to n, while they must move move to the right thereafter, 
returning to their starting points at time 2n. 

To count the configurations, label the particles by their initial location of the one-dimensional 
line I = 1, . . . ,p. At each time step one walker will move one lattice site to the left (L) or right 
(R), subject to the rule that no two particles can occupy the same lattice site (in the counting 
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Figure 1: A particular configuration of 3 random turns walkers performing 8 steps in the sequence 
L 4 i? 4 . The walk can conveniently be represented diagramatically as done at right according to the rules 
specified in the text 



problem we take w\ = The constraint that the particles return to their initial positions 

after 2n steps requires that for each walker the number of steps to the left equals to the number 
of steps to the right after 2n steps, and that in total there are n steps L and n steps R. This 
latter fact allows the walks to be partitioned according to the ordering of L's and i?'s, of which 
there are y™j possibilities. The simplest of these is L n R n , which means the first n steps are 
all to the left, while the final n steps are all to the right. We now pose the question: for a 
given combination of {L n , R n } corresponding to a particular ordering of L, R steps, how many 
distinct configurations of the p random turns walkers are there? 

Consider the particular ordering L n R n . We will represent each configuration diagramatically 
by placing a square marked t in column k (k = 1, . . . ,p) if walker k moves at time step t. This 
square must be placed above any other squares in the column if the step is L, and below any 
other squares in the column if the step is R. If there are no other squares in the column (i.e. the 
walker k is making its first move) then the square is placed immediately above the axis if the 
step is L, or immediately below the axis if the step is R. An example of such a diagramatical 
representation is given in Figure |]. 

The problem now is to determine the number of distinct diagrams. Consider the region 
of the diagram above the axis. Here each column, of which there are no more than p, must 
be weakly decreasing in height. Each of the n squares must be labelled by a different integer 
l,...,n with the numbers strictly increasing along rows and up columns. We recognize such 
a diagram as equivalent to a standard Young tableau (see e.g. [11|; the Young tableau results 
be interchanging rows and columns). Thus there is a bijection between the number of walks 
following the sequence L n and standard Young tableaux with entries 1 , . . . , n and no more than 
p rows. We remark that this is not the first time a bijection between Young tableaux and vicious 
walker problems has been observed: in |jl2[ a bijection between semi-standard tableaux and the 
configurations in a special case of the lock-step random walker model || || was identified. 

Regarding the last n walks, because each walker must return to its starting position after 
the 2n steps, the shape of the diagram below the x-axis must be the mirror image of the shape 
of the diagram below the x-axis. Furthermore, by pushing each column k = 2, . . . ,p of this 
downwards so that they all align with the first entry of the first column, and then relabelling 
each box by t i— » In + 1 — t we see that another standard Young tableau results. (This procedure 
is equivalent to constructing this diagram starting with the final step of the walk and working 
backwards.) These transformations are indicated in Figure [2]. 
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Figure 2: Transforming the diagram below the horizontal axis in Figure [j] into equivalent Young tableau 
form 



Thus configurations of the p random turns walkers under consideration are in bijective corre- 
spondence with pairs of standard tableaux of n-boxes constrained so that there are no more than 
p rows (recall the role of rows and columns in Young tableaux and our diagrams is reversed). 
But such pairs of standard Young tableaux are well known (see e.g. fliT| ) to be in bijective corre- 
spondence with permutations of {1, . . . , n} such that the length of each increasing subsequence 
is less than or equal to p. Hence we have enumerated the number of walks in terms of such 
permutations. 

Proposition 1 Consider p random turns walkers, initially equally spaced one unit apart and 
returning to their initial position after 2n steps. Suppose the walkers make their steps in the 
sequence L n R n . The total number of distinct configurations equals the number of permutations 
of {1, . . . ,n} such that the length of the maximum increasing subsequence is less than or equal 
to p. 

Consider now another sequence of n L's and n i?'s. This sequence can be transformed into 
the sequence L n R n by elementary transpositions Si which interchange the ith and (i + l)th 
members of the sequence, assumed to be R and L respectively. Likewise, we can define the 
corresponding action of Si on a diagram (or equivalently the lattice paths) and so obtain a 
bijection between distinct configurations with walks following the sequence L n R n , and distinct 
configurations with walks following some combination of the sequence L n R n . We assume step i 
is opposite in direction to step i + 1, and defined the action of Sj to first interchange the position 
of boxes i and i + 1. At the level of the lattice paths this has the action depicted in Figure |3[ 

If the two boxes are in the same column, it may happen that the resulting lattice path is 
inadmissable, in that the new (local) configuration intersects with an existing path (note that 
this cannot happen in the first two cases of Figure In such a circumstance we move the 
configuration to the left (right) in the second last (last) situation of Figure |3| until a permissable 
configuration is obtained (after so moving the left-right pairs two vertical lines, corresponding 
to a stationary walker, take its place). This is illustrated in Figure ||. Notice that in all cases 
sf = 1, so the procedure is invertable. Furthermore, the braid relations SjSj+iSj = Sj+iSjSj+i are 
satisfied so the correspondence is independent of the order of application of the transpositions. 

Let's now describe this aspect of the definition of the action of Si at the level of the diagram. 
For this purpose we must specify an admissable diagram. From the rules of the random turns 
model we see that an admissable diagram is constructed by putting boxes in the columns 1 to p 
subject to the conditions that box i must go above (below) the axis if the step is L (R), and that 
for each column the number of boxes above the axis minus the number below must be greater 
than or equal to the same quantity for the column on the right. Some examples of such diagrams 
and the corresponding combinations of {L n , R n } are given in Figure ||. So we want to describe 
the situation in which boxes i and i + 1 are interchanged but an inadmissable diagram results. 
According to Figure |] both boxes should be shifted to the left (right) if the new position of box 
i is below (above) the axis until an admissable diagram results. Also we adopt the convention 
that numbers in each column should be increasing above and below the axis, so we rearrange 
the boxes within a column above and below the axis accordingly. 
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Figure 3: In the first two cases the distinct walkers correspond to the columns of boxes i and i + 1, 
while in the last two cases boxes i and i + 1 are in the same column, and the resulting configuration 
assumed admissable 



Figure 4: Here both box i and i + 1 are in the same column and the corresponding walker configuations 
going from step i to step i + 1 are as indicated. Because the rule of Figure || would lead to inadmissable 
configurations in each case, Si acts by propagating the new configuration to the left and right in the two 
cases respectively until an admissable configuration is obtained. 
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Figure 5: An example of the action of the elementary transposition operator on some permiss- 
able diagrams in the case that both boxes are from the same column with the situation of Figure 
applying. The corresponding sequences of L's and i?'s are also given. 
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The action of the elementary transpositions thus described give a bijection between lattice 
paths with steps in the sequence L n R n , and lattice paths with steps in a sequence of any 
particular combination of {L n ,R n }. Thus by making use of Proposition [l] we can solve the 
enumeration problem for all such walks. 

Proposition 2 The result of Proposition [J for the number of walks in the sequence L n R n also 
applies for any combination of {L n ,R n }. 

It is of interest to note that there are multiple integral formulas for both the number of 
random permutations of {1, . . . ,n} with at most p increasing subsequences, and the total number 
of random turns paths with p walkers starting at sites l[, . . . ,l' p and finishing at sites li,...,l p 
in 2n steps. Let us denote these numbers by f np and Zmtfx, ■ ■ ■ , lL\ hi ■ ■ ■ > lp) respectively. Then 



we have [19] 



^tS(^/>--1>(£ 2 ^) ! " n .v-*? (i) 



j=l 1<"</3<P 



and U 

Z 2n (l[,...,l';l 1 ,...,l p ) = -— (* Mi- r deJY / 2co S e j ) 2n det\e-^- l 'e )e ' . 

(2) 

Let us specialize (||) to the situation of Propositions [j] and [2] by choosing lj = lj = j 
U = 1, • • • ,P)- Now 

det[e-^~^U =1 _ p = f[ e"*- 1 *' det^- 1 ^]^,...,, 

3=1 

= JJe-*^" 1 ^ J] (e^-e^), 

j=l 1<"</3<P 

where the final equality follows from the Vandermonde formula. Since the product of differ- 
ences is antisymmetric the non-symmetric factor in the integrand nf=i e - ^-? -1 ^ can be anti- 
symmetrized, giving another Vandermonde product, provided we divide by pi. Thus we have 

Z2n({l'j = j}j=l,...,p', {h = k}k=i,..., p ) 

1 1 r d6 X --- r d0j^2cos%) 2n n ^^-^f. ^) 



T j=l 1<«</3<P 



p\ (2tt)p 
Comparing (||) and (|3|) gives 

Z2n({l'j = j}j=i,..., P ; {h = k}k=i,..., p ) = ( 2n )/n P (4) 



which is of course also an immediate corollary of Proposition |T| and |2l However once having 
deduced (H), the formula (|^) for f np follows as a special case of (|2|). 

We know from the derivation of Proposition || that there is a bijection between configurations 
of p random turns walkers performing 2n steps in the sequence L n R n before returning to their 
initial positions of all one unit apart, and pairs of standard tableaux each of the same shape 
and consisting of n boxes. Furthermore, in the bijection the length of row j corresponds to the 
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maximum displacement of walker j to the left of its starting point (this occurs at step n). Thus 
if we choose p > n this same bijection holds but with no restriction on the number of rows. In 
such a situation the asymptotics of the row lengths are known precisely jjj [18|, §,14]. We can 
therefore give these results an interpretation in the random walker setting. 

Proposition 3 Denote the positions of the walkers on the one- dimensional line in the above 
asymmetric random turns model by L where initially walker j is at position lj = j. Define the 
scaled displacements by 

^W/3(_|__ 2 ) (5 ) 
and the corresponding scaled k -point distribution function by 



Ik) := Jim (-y<-) k Pk\h, . . . , f fc ) 



Pk{ 

n-»oo \ n * 

where p) denotes the k-point distribution for the walker problem in the finite system. Then 
from the results of ]Tj| , \L% ] for the tableau problem we have 

Pk{h, ...Jk) = det K(l a , lp) (6) 

L Ja,p=l,...,re 

where 

Ai(s)Ai'(y)-Ai'(s)Ai(i/) 
K(x,y) := — . 

The distribution (JgT) is precisely the scaled distribution of the eigenvalues at the edge of the 



spectrum for GUE random matrices [1C, ^0[. Perhaps more relevantly to the random walker 
problem, (|5|) coincides with the scaled distribution for free fermions on a line confined by a 
one-body harmonic potential, at the edge of the support of the density. The relevance is that 
there is a well known relationship between continuous models of non-intersecting walkers and 
free fermions (see e.g. ||). 

Regarding some physical features of Proposition |j| note from that the average displace- 
ment is fx = In 1 ! 2 (for a recent independent proof of this result see |l~3|l), with standard deviation 
proportional to (An) 1 / 6 = /i x , \ = 1/3- As emphasized in the exponent x = 1/3 is typical 
of two-dimensional growth models (it can be derived from the one-dimensional Burgers equation 
describing such processes pl|) . On this point we recall that vicious walker paths fixed at the 
endpoints as in Proposition^ form the well-known (see e.g. |l6|]) terrace-step-kink model of a 
crystal surface. 
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